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Abstract: We study the ground state energy of integrable 1 + 1 quantum field theories 
with boundaries (the genuine Casimir effect) . In the scalar case, this is done by introducing 
a new, "R-channel TBA", where the boundary is represented by a boundary state, and 
the thermodynamics involves evaluating scalar products of boundary states with all the 
states of the theory. In the non-scalar, sine-Gordon case, this is done by generalizing 
the method of Destri and De Vega. The two approaches are compared. Miscellaneous 
other results are obtained, in particular formulas for the overall normalization and scalar 
products of boundary states, exact partition functions for the critical Ising model in a 
boundary magnetic field, and also results for the energy, excited states and boundary 
S-matrix of 0{n) and minimal models. 
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1. Introduction 



The main purpose of tiiis paper is to study the ground state energy of 1+1 integrable 
relativistic quantum field theories with boundaries. This involves several questions. One 
is the energy associated with a boundary for an infinite system, the other is the way the 
energy of the theory on an interval varies with its length - the "genuine" Casimir effect. 

Some of these properties are most easily studied directly in a continuum formalism, 
using a scattering description for either a massive or massless bulk theory, the information 
on the boundary being encoded into a reflection matrix. Other properties are most easily 
studied using a lattice regularization. 

In computing energies, we shall use detailed information about the boundary states, 
and this work also presents new results on that question. For instance, the "R channel" of 
the Thermodynamical Bethe Ansatz (TBA) approach provides expressions for the scalar 
products of boundary states. 

For a Quantum Field Theory defined on a torus, the standard way to compute its 
ground state energy is through the Thermodynamic Bethe Ansatz [|]. If the theory is 
defined on a circle of circumference R, one switches to a modular transformed point of 
view where now the theory is defined on a circle of very long circumference L and at 
temperature 1/R. The free energy F of the theory in the "i? channel" can be computed 
using TBA, and it is simply related to the ground state energy of the theory in the 

"L channel" by F = -TLE^{R). 

Consider now a quantum field theory on a cylinder of finite length R and circumference 
L, with some boundary conditions (a, h) at the ends of the cylinder. There are two possible 
ways of viewing the partition function: in the 'R-channel' where time evolution is in the 
R direction, and in the 'L-channel' where time evolution is in the L direction. (See figure 
1). This problem is more difficult because the theory looks rather different in the two 
channels. To use the same argument and do a TBA in the R channel, one now needs to 
control the effect of boundary states on the thermodynamics. This is discussed in section 
2 for the simple case of the Ising model to clarify a few delicate points, and in section 3 
for the general case of a scalar theory. The TBA in the L channel gives rise to a different 
kind of information, related to the normalization of boundary states. This is discussed in 
section 4. Finally, all these results can be generalized to the case where the theory is still 
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Figure 1. Cylindrical geometry in the R and L channels. 



conformal in the bulk but has a boundary interaction that breaks the conformal invariance 
(section 5). 

In the second part of this paper we address the same questions for a non scalar theory, 
namely the sine-Gordon model with Dirichlet boundary conditions. This case is technically 
more difficult. Instead of rederiving a TBA in the R channel, we have generalized the 
beautiful approach of Destri and de Vega for a system with boundaries. This requires 
working for a while with the lattice theory, here chosen to be the XXZ model with boundary 
magnetic fields. 

In section 6 we recall the various elements of the model, and carry out the TBA in 
the L channel. As a result, we are able to obtain the boundary ground state energy of the 
sine-Gordon model with Dirichlet boundary conditions. In section 7 we write the ground 
state energy. In section 8 we use this result for minimal models, proposing in particular 
boundary conditions that correspond to excited states of the theory. Section 9 collects 
conclusions. In particular we show that the DDV equations are actually essentially similar 
to the equations for a scalar theory derived in section 3. 
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Part I: Scalar theories 



2. Ground state energy of the Ising model with boundaries 

In the bulk, the Ising model can be described in terms of massive fermionic operators 
A{e) and A^O) with the usual anti-commutation relations, {A{9), {9')} = 2'k5{9 - 9'), 
where 9 is the rapidity variable. The mass m of the fermionic field is proportional to 
Tc — T. The S'-matrix in the bulk is simply —1. In this section, we are interested in 
the explicit computation of the ground state energy of this model for different boundary 
conditions. Let us consider then the Ising model in the low temperature phase (T < Tc) 
in the geometry of figure 2, i.e. a long strip of horizontal length L and vertical width R. 
At the extremities of the horizontal axis, the order parameter of the model is subjected to 
periodic boundary conditions while in the vertical direction it satisfies boundary conditions 
of type a and b. Thus, we are effectively considering the partition function on a cylinder 
of circumference L and length R, with different boundary conditions on each end of the 
cylinder. (See figure 1.) As we discuss in more detail below, such boundary conditions are 
expressed in terms of boundary states | Ba > and \ Bf, >. 
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Figure 2. Strip geometry. 
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Depending of our quantization scheme, we have two possible ways to compute the 
partition function of the system. The first possibihty consists in choosing as direction of 
time the horizontal axis and therefore the partition function will be expressed as 



Tre-^^-\ (2.1) 



where Hab is the Hamiltonian relative to the system with boundary conditions (a, 6). In 
the second method the time evolution takes place along the vertical axis and therefore the 
partition function is given by the matrix element of the time evolution operator between 
the boundary states, i.e. 

Zab =< Ba I e-^^ I > , (2.2) 

where now H is the Hamiltonian of the bulk system. The ground state energy Eab{R) is, 
by definition, the leading term arising in the large L limit of the first expression, eq. ( |2.1|) , 
i.e. 

Zab - e-^^-(^) . (2.3) 

However, in view of the equivalence of the two quantization schemes, we can compute this 
quantity by looking at the large L limit of the second expression, eq. ( |2T^ ) . 

Since eq. ( pT^ ) employs the boundary states of the model, let us shortly recall their 
basic properties (for more detail see the original reference 0). The Ising spins placed on 
the boundaries can be subjected to three possible boundary conditions, namely: (i) they 
can be frozen to one of the two possible fixed values ± (fixed boundary condition); (ii) they 
can be completely free to fiuctuate (free boundary condition); (iii) they can be coupled to 
a boundary magnetic field h (magnetic boundary condition). In the QFT description of 
the model, each of the above microscopic configurations corresponds to a boundary state 
I S >, which for infinite length L, reads 



B >= gexp 



^Kie)A\-e)A\e) 

Ztv 



> . (2.4) 



Here g is an overall normalization which we shall discuss later. Since we are for the moment 
interested in the large L limit, we simply set this factor equal to unity. For simplicity, 
we ignore possible additional contributions to the boundary state from zero momentum 
particles. From the point of view of QFT, the boundary state can be therefore regarded as 
a particular state of the Hilbert space of the bulk theory, made of a superposition of pairs 
of particles of equal and opposite momentum ("Cooper pairs"). All information relative 



4 



to a particular boundary condition is encoded into the function K(9) which can be seen 
as the elementary amplitude to create a virtual pair of particles. 

For the Ising model at low temperature, the vacuum state | > can be one of the 
two possible vacua of the model in this phase and the explicit expressions of the amplitude 
K{9) for the three cases above considered are given by 



for the free boundary conditions, 



i^/ree = -«COth - , (2.5) 



-ft^/ixed = « tanh - , (2.6) 



for the fixed boundary conditions, and 

, 6 k + cosh 9 

for magnetic boundary conditions, where k = 1 — It is evident that by varying /i, we 
can interpolate between the free boundary condition {h = 0) and the fixed one [h oo). 
After this brief discussion on the boundary states, let us come back now to the eval- 



uation of eq. (|2.2| ) in the limit L ^ oo. The first thing to consider is the action of the 
time evolution operator e~^^ on the boundary state \ Bb >. Since H is the Hamiltonian 
of the bulk theory, it is diagonalized on the basis of the multiparticle asymptotic states 
I 9i,...,9n >, with eigenvalues given by m^^cosh6'i. Expanding the boundary state 
\ Bb > in terms of multiparticle states by using its definition ( |2.4|) , it is easy to see that 
the action of the time evolution operator on | > simply reduces to a redefinition of the 
amplitude Kh{9) of this state, namely 

Kb{9) ^ Kb{9, R) = Kb{9)e-^'^^^°^^^ . (2.8) 

Consequently, from an abstract point of view, the evaluation of the partition function Zab 
consists in computing the scalar product of two boundary states, one of them R dependent, 
i.e. Zab =< Ba I Bb{R) >. To determine such a quantity, it is convenient to expand both 
boundary states in their multiparticle components. For the integrability of the theory, 
transition amplitudes between states with different number of particles are not allowed 
and therefore Zab can be expressed as 



7V=0 
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where 



< I A{-9[)A{e[) . . . A{~9'j^)A{9'^)A\9n)A\~9n) • • • A^{9i)A\-9i) \ > 

(2.10) 

The first term of this series is given by Z^^^ =< | >= 1. The higher terms can be 
computed by using the Wick theorem. There is, however, a subtlety which already arises 
in the computation of the second term 

2^= Z ^K^m,{9',R){2nr6{9-9')6{9-9') = 
Jo 271^0 2n 

roc iQ_ ■ l^-J-J-J 

= 5(0) / ^i^„(^)i^,(^)e-2-^'=°^^^ 
Jo 27r 

The term 5(0) signals the divergence due to the infinite volume limit on which the boundary 
states are defined. To correctly extract the L dependence of this expression, we need to 
redefine the normalization of the multiparticle states, i.e. to provide a regularization of 
the 5 function. If we were working in the momentum space, a possible regularization of 
the 5 function is given by 

{p) = ^ dx e'P^ = — sin ^ . (2.12) 
2n J_L Tip 2 

Using this definition, we can make sense of the square of the 5 function and in the limit 
L — > oo we have {5l{p — q))"^ — ^ 'i^^iP ~ Q)- However, working in the rapidity variable, 
the final expression of the square of the d function entering ( 2.11 ) acquires an extra term 



m cosh 9 due to the Jacobian of the transformation from the momentum to rapidity variable 
and the final result is to substitute 

[6(9 - 9')f ^ — cosh 9 5{9 - 9') . (2.13) 
27r 

This regularization allows us to keep track of the L dependence of our expressions and 
therefore to extract the ground state energy E^^^. With this substitution, Z^^^j is now 
given by 

poo in 

Z^^^=mL ^cosh9Ka{0)Ki,{9)e-^'^^''°'^'^ . (2.14) 
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For later convenience, it is useful to define the following quantities 

r ^coshe[Ka{0)Kh{e)e-^'^^^°'^^]^ . (2.15) 
Jo 27r 

With the new normalization of the states, the higher terms Z^^^^ in the series ( p.9|) are 
polynomials in the variable (mL) of order A^, as explicitly shown by the first representatives 



(4) _ {mLf 2 (wL) 

(6) {mLf 3 [mLf (mL) 
^ab = 3, h 2 — ^ ^ ' (2.16) 

(mL) 



(8) _ (mL) 4 (^n^T^T ^ ^""^^ 
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and so on. The resulting sum organizes into an exponential series 

Z^, = l + {mL)f{R) + i!^(/(i?))2 + . . . i!I^(/(i?))- + . . . = exp[mL/(i?)] (2.17) 
where 

f{R) = = / 1^ cosh ^ log (l + ¥„(^)i^,(^)e-^-«^°^^^) . (2.18) 

N=l 

Therefore, the ground state energy for boundary conditions (a, b) is given by 

/oo 
ci^cosh^log(l+¥„(6')i^fe(6')e-2'^^'=°^^^) . (2.19) 
-00 

The above result can also be easily derived as a special case of the general formula 

< o|e("'^^")e("^'^°^^)|0 >= det(l + A^M), (2.20) 

where (a. Ma) = J2n,m^riMnmam and {an,aj„} = {an,aj„} = 5n,m, {an,aj„} = 0. Spe- 
cializing to our situation, the matrices M^N are kernels: M{9,9') = 6{9 — 9')K{9), 
N{9,9') = d{9 - 0')i^(^)e-2'^^'=°"^^ and the determinant in ( ^M) is a Fredholm de- 
terminant. Using logdet(l -|- NM) = trlog(l -|- NM), and regulating the squares of delta 
functions as in ( |2.13| ), one obtains ( [2.191 ). 



As a check of this expression, we can compare its ultraviolet limit R ^ with the 
result directly obtained by using CFT H]. In the conformal limit, general expressions for 



the partition function in the cyhndrical geometry with boundaries were derived by Cardy 
1^. In the R-channel, one has 

Zab = Yl<bXr{q), (2.21) 

i 

where q = exp(— 7rL/i2), Xi are the characters, XiiQ) = q~^/'^^Tr , and denotes how 
many times the representation [i] appears with the boundary conditions (a, b). Thus, as L 
tends to infinity, one has 

= I [lT - ^) , (2.22) 



R\ 24 

(0) 



where Lq is lowest Lq eigenvalue of the states in the representations [i] for which n\^^ ^ 0. 



Taking the limit — in (|2.19|) , the ground state energy assumes the scaling form 

1 

E'abiR) ^ delog [1 + Xab{oo)e-^] , (2.23) 

where 

Xab{e) = Ka{e)Kt{e). (2.24) 

Since Xfreejreeioo) = X fixed Jixedioo) = 1 , 

E free, free -^/jxedj/ixed ^ 48-R (2.25) 

With those boundary conditions, the ground state energy is dictated by the conformal 
family of the identity operator. For the mixed boundary conditions, we have instead 
Xfree,fixed{oo) = -1 and therefore 

E%ee, fixed ~ • (2.26) 

In this case, the ground state energy is ruled by the conformal family of the magnetization 
operator, of conformal dimension j^. These results are in agreement with those obtained 
by Cardy M. 



3. The ground state energy of a scalar theory with boundaries and TBA in 
the R-channel 

Reproducing the above computation in the case of a scalar interacting theory is more 
difficult. This is because the particles scatter non trivially, so care must be exercised in 
computing scalar products. Also, the allowed rapidities of the physical states are non 
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trivial solutions of Bethe type equations. The simplest way to proceed is to use a "thermal 
approach", generalizing the standard Thermodynamic Bethe Ansatz (TBA). 

The TBA equations in a torus geometry with no boundary were studied in . In the 
cylindrical geometry considered here with boundary interactions, the two possible ways 
of viewing the partition function which correspond to the two possible channels for time 
evolution are rather asymmetrical, thus one expects the TBA equations for each channel 
to take different forms. 

In the analysis of the previous section, we examined the large L behaviour of the 
partition function using the boundary states. This corresponds to viewing the partition 
function in the 'i2-channel', where time evolution is in the direction of R. As we now 
describe, one can obtain the TBA equations in this channel in the limit of large L. 

In the situation with non-trivial S-matrix, one may attempt to evaluate (|2.2| ) using 
the form of the boundary states ( p.4| ) and the Faddeev-Zamolodchikov algebra: 



A\e)A\e') = s{e - e')A\e')A^{e) , 

A{e)A{6') = S{6 - 6')A{6')A{6) , (3.1) 
A{6)A\6') = S{e - 6')A\6')A{e) + 5{6 - 6') . 

As in the Ising case discussed in the previous section, one finds that the result is ill-defined 



due to the fact that the formula ( |2.4| ) for boundary states is meant in the limit L — oo, 
and this leads to squares of 5-functions. In the Ising case, the regularization ( p.l3|) is 
adequate since the particles are free and the quantization condition on the momenta does 
not involve the S-matrix. More generally, one can view the following TBA derivation as a 
way of regularizing the inner product ( |2.2|) for very large but finite L. 

For simplicity, consider a model with a single particle excitation. Our starting point 
is eq. ( p.2| ) which can be written as 

Z^^J2 <Ba\a><a\B,> ^_^,^^ (3.2) 

< a\a > 

a ' 

where \Ba > and \Bi, > are boundary states and the sum over a stands for a formal sum 
over all the states of the theory. Due to the form of the boundary states, the states \a > 
that contribute to the sum are of the form 

(3.3) 

= A\eN)A^{-eN) . ..A^{9i)A\-9i)\0 >, 
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where 6^ > On-i > ■ ■ ■ > 9i > (inequalities are strict, since due to the condition 
S{0) = —1, the Bethe wave function for particles on a circle of radius L vanishes when two 
rapidites are equal). Using the Faddeev-Zamolodchikov algebra, one finds that 



N 



<Ba\2N>= (5(0))^n^«(^^)' (3.4) 



i=l 



where the formal expression (5(0))"^ results from terms 5{6i — 6i). It is simple to show 
(U) by introducing pair creation operators B{6) = A{-0)A{6), B\e) = A^{0)A^{-0), 
satisfying the algebra 

B{9)B{e') = B{e')B{e) , B\e)B\e') = B\e')B\e) , 
B{e)B\e') = B\e')B{e) + 5^{e - e') (3.5) 
+ 5{e - e')s{e' - e)A\-e')A{-e) + 5{e - e')s{e - e')A\e')A{e) . 

One also finds < 2N\2N >= {S{0))^^, thus, 

1=1 

Let us introduce a density per unit length of pairs of particles P{0), so that in ( |3.2|) , 
Ea = '^^^i 2m cosh 6 i is replaced with 2mL J d6 cosh.{6)P{6). Due to the factorization 
( ^^) , the partition function reads 

Z ^ j[dP]expS^L [log(K a{0)Kb{9))-2Rm cosh e]P{e)d9 + S{[P])Y (3-7) 

where S[P] is the entropy of the particle configuration described by the distribution P{6). 
From this equation, we see that the contribution of the boundaries has a natural interpre- 
tation as (rapidity dependent) chemical potentials. 

The density P{0) is constrained by the quantization condition. For the states |2A^ >, 
this reads 

e^mLsir.\.{e,)g^26i)WS{6i-ej)S{6,+6j) = 1. (3.8) 

Introducing densities of particles and holes as usual such that L[P{6) + P^{6)]d6 is the 
number of allowed rapidities between 6 and 6 + dO, one has the quantization conditions 

poo 

2ti {p{e) + p'^ie)) = m cosh e-27t [^{6 - e') + $(6* + e')] p{e')de' , (3.9) 
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where $ is defined by 



1 d 

2inde 



iiiS{e). 



(3.10) 



In terms of P, P^, the entropy takes the usual form: 

/"OO 

S{[P]) = L de [{P + P'') \n{P + P^)-P In P-P^ In P''] . (3.11) 
Jo 



Finally, the TBA equations arise as a saddle point evaluation of (|3.7| ), i.e. the minimization 



of the free energy, subject to the constraint (|3.9|). Using standard manipulations, one 



mL 



-e L: 






/"OO 






/ cosh 9 In 


de 


/o 







where e satisfies 



e(^) = 2Pmcosh^ - log{KaKb) + * log(l + e" 



(3.12) 



(3.13) 



In the previous equation, $s = $(6* — 9') + ^{9 + 9') and * denotes the convolution. The 
combination Ka{9)Kb{9) is an even function of 6*, since K{9) = K{—9). Thus, we can 
extend the domain of definition of e(^) to the whole real axis by t{—9) = e{9) and let the 
integrals run from — oo to +oo. The equations for the ground state energy may be recast 
into the following: 



cosh 9 log 



1 + Xab{9)e 



d9. 



e{9) =2i?m cosh 9+ ^{9 -9') log 1 + Xab{0')e 



(3.14) 



d9', 



where 



Xab{9)=Ka{9)Kb{9). 

One may now see explicitly how the TBA regularized the inner product ( |2.2| ). A 
multiparticle expansion of the TBA equations ( |3.12| ) ( pT3| ) can be obtained by solving 
them iteratively as follows: 



e^+i = eo + $. *log(l + e '") 
starting from eo = 2mi?cosh6' — log{KaKb). One obtains 

Z,, = l + —{D^+D2 + . ..) + -( — ) Dl + 



(3.15) 



(3.16) 
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where 

/oo 
d6 cosh6e-^"'^''°'^^KaiO)Kb{e) , 
-OO 

-| poo -| poo 

D2 = / decOshee-'^"'^'='''^^Kl{e)Kl{e) + — ^Q^Q'^-^mRicoshe+coshB') 

($(^ + 0') + ^0 _ 0'))Kaie')Kb{e')Ka{e)Kb{e). 

(3.17) 

On the other hand, evaluating ( p.2| ) using the Faddeev-Zamolodchikov algebra, one obtains 



+ dede's^{o)Ka{e)Kt{e)Ka{o')Kt{e')e-^'^^^''°''^^+''''^'^'^ + .... 



(3.18) 



Comparing (|3.16| ) with (|3.18| ) , one sees that the single 5-functions in ( p.l8|) were regulated 
as in ( |2.13| ), however the 5^(0) is regulated in a more complicated fashion involving the 
S-matrix: 

(J- \ 2 J- 
^) cosh ^ cosh 6*' + ^ cosh ^($(^ + ^') + ^(^-^'))- (3-19) 

As an explicit example of our discussion, we will consider the 991,3 massive deformation 
of the minimal non-unitarity conformal model A^3,5, in the R-channel. The central charge 
is c = —3/5 and there are four primary fields of chiral dimensions (0, — ^,;|,|). The 
massive phase of this deformation has only one massive particle with the two-particle elastic 
S'-matrix given by S{9) = — ztanh ^ {O — i^) @. For the boundary reflection amplitudes, 
solutions of the equations 



n{e)n{-e) = i, 

we have two minimal solutions, given by 

sinh(f-f) 



(3.20) 



ni{e) = -itanh 
7^2(^) = zcoth 



2 4;sinhf + f ' 

. y 3.21 
i7r\ sinh(f - f ) 



2 A ) sinh (f + f ) ' 

and from those, we can compute K[9) = 71(^^ — 9). We now want to extract the ultravi- 
olet limit of the ground state energy for different boundary conditions associated to these 
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amplitudes. The kernel entering the TBA equation is given by $(0) = 1/ cosh^. As in 
analogous TBA computation, the ultraviolet limit of the ground state energy is ruled by 
the "kink" solution of the pseudoenergy e(^), i.e. the universal function which flattens in 
the central region — log(2/r) <^ 9 <ti log(2/r) to the constant value e^^^, solution in this 
case of the algebraic equation 



.0 



-^log l + Xab{oo)e 



(3.22) 



In fact, in the ultraviolet limit mR — ^ 0, we have 

1 r , ( eXab(oo)e-' 



AtiR 



de 



1 + Aa6(oo)e" 



+ log [1 + A„b(oo)e-^] 



(3.23) 



Choosing boundary conditions of the type (1,1) or (2,2) on both sides of the strip, we 
have Aii(cx3) = A22(oo) = 1 and for the constant value of the pseudoenergy we have 



-11 



"22 



log 



l + v^' 



(3.24) 



and final result 



^?i 



^22 



TT 



(3.25) 



On the other hand, choosing boundary conditions of the type (1, 2), we have Ai2(oo) = — 1. 
In this situation we have 

, ,'-l + v^' 



"12 



and correspondingly 



log 



eUR) 



TT 



(3.26) 



(3.27) 



In terms of CFT, in the first two cases the boundary conditions select as dominant con- 
tribution in the partition function the one of the conformal chiral field a of anomalous 
dimension A = — ^ whereas in the last case, the leading contribution comes from the 
identity family. 

In closing this section, observe that by taking the limit i? — in the above equations 
( P.12|)(|3T3D , we obtain the leading behaviour for the scalar product of boundary states 



In < Ba\Bb >^ 
e{9) = 



mL 



4tv 



cosh 9 In 



1 + Xabe 



^{9-9')ln 1 + Xab{9)e 



<0) 



d9, 
d9. 



(3.28) 
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4. Normalizations and TBA in the L-channel 

We can also discuss in more detail the overall normalization alluded to in ( |2.4| ). There, 
|0 > is the ground state of the massive theory, which in the renormalized theory is simply 
a state without any excitation. Although the correspondence between particle states of 
the massive theory and conformal states in the UV limit is not totally known, the ground 
states usually match, that is in the UV limit, |0 > in (|2.4| ) goes to the conformal ground 
state |0 >c, while all the states with particles go to combination of fields and descendants 
with higher conformal weights. Therefore, if the boundary state ( |2.4| ) corresponds to some 
state \B >c in the conformal theory, we expect 

g =< 0\B{m) >= c < 0\B >c for any m. (4.1) 

Let us further investigate this identity. In the "L-channel", time evolution is along the 
circumference of the cylinder of length L, and the Hilbert space is constructed along the 
finite segment of length R. In the limit of large R, TBA equations for the partition function 
were obtained in . In this channel, the particles are confined to the segment of length R 
and refiect off the boundary with refiection amplitude TZa{9) at one boundary, and TZb{d) 
at the other boundary. The quantization condition on the momenta involves in this case 
the boundary refiection amplitudes: 



g2jmi? sinh(6i) 



7^a(^^)7^6(^^) = l, (4.2) 



where 9i > 0. Introduce the density of particles and holes as usual, so 2R[P{9) + P'^{9)]d9 
is the number of allowed rapidities between 9 and 9 + d9. It is also convenient to let 
rapidities run from — oo to oo, defining the densities for negative rapidities by parity, so 
(U) reads 

27T{P{9)+P^i9))=mcosh9-2n [ ^{9 - 9')P{9') + , (4.3) 

where $ is as before, and 

Qab{9) = i A log(7^„7^,) - -^\ogS{29) - 27:5{9) , (4.4) 

I CLu I CLu 

and we used the fact that, by unitarity, Oab is an even function. In the latter expression, 
the 5 term arises because vanishing rapidities, although possible solutions of (|4.2|) , are not 
acceptable since the corresponding wave function vanishes exactly. 
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From here the standard TBA procedure leads to the foUowing result: 



1 

log Zab / [2mi?cosh6' + QabiO)] log 

/oo 
^{6- 6') log 
-oo 



l + e" 



1 + e- 



d6, 

(4.5) 



de' 



The © term in ( [4.5| ) is i? independent and may be considered as a boundary free energy 
term. As discussed in the TBA actually reproduces the boundary free energy up to 
an additive constant. This is because, when considering contributions of order one, one 
should take care of various corrections like corrections to the Stirling formula in evaluating 
the entropy, or loop corrections to the saddle point, which we have discarded here. Also 
the correspondence between the entropy of the field theory and the one computed using 
particles might involve some constant, for instance if the particles describe kinks and an 
overall choice of configuration is left out. 

Now we can interpret this boundary free energy easily. Indeed in the large R limit, 
the next to leading behaviour of the partition function (|2.2|) is 



Z„b ^< S„|0 X 0|S, > e-«^o(^), (4.6) 

up to exponentially small terms, where Eq is the ground state energy of periodic Hamil- 
tonian. Therefore we have 



oo 



In < SJO >< OlSb >= -3- / e„6(^)ln l + e-'(^) d^ + cst. (4.7) 

In the IR limit given by L ^ oo, e goes to infinity and the first term on the rhs of ( [4.7| ) 
vanishes. In the UV limit, expressed by L 0, on the contrary e goes to a constant. The 
statement (|4.1|) translates therefore into the condition 

Q{e)de = Q. (4.8) 



This condition is indeed true, as can be easily seen by using (|4.4| ) and the cross unitarity 
relation 

n{i^-e)=s{2e)n(i^ + e), (4.9) 

together with analyticity. 

As we will see in the next section, for finite L, g can be a non-trivial function of L. 
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5. Massless Theory in the Bulk 

An interesting situation is when the bulk theory is massless, but the conformal invari- 
ance is broken by the boundary interactions. The previous formalism carries over easily 
to that case. This time the bulk theory is described by massless particles which are left or 
right moving with trivial left-right scattering. Parametrize the right (left) moving particles 
by rapidities so that E = p = (resp. E = —p = ^e~^) (observe that opposite mo- 
menta still correspond to opposite rapidities) . Then the left- left and right- right scattering 
are described by an identical S matrix in terms of rapidities, and one can introduce as 
before 



m 



1 d 



In Sll{9). 



2m dO 

Now all rapidities in [— cxo, oo] are allowed, so the quantization condition reads 



m 



2'K{p{e) + p\e)) = -e^ 



2n / ^{e - e')p{e')de' , 



(5.1) 



(5.2) 



and for the entropy ( p.ll|) the integration runs from — oo to oo. The energy reads E = 
YliLi "le^' . To have a breaking of the conformal invariance and a corresponding flow, 
one needs to put an energy scale at the boundary, so the boundary S matrix depends on 
some extra parameter (e.g. a magnetic field), corresponding to an amplitude Ka{9 — Ob^)- 
The boundary state is then a combination of massless states involving pairs of left and 
right moving particles with opposite momenta (hence opposite rapidites using the above 
parametrization) 



\Ba >= fi'exp 



d9 
2^ 



Ka{e-eBjAi{-e)Ai{e) 



10 > . 



(5.3) 



Then 
to 



still holds. In the "R-channel" , the maximization of the free energy leads now 



E': 



m 
4tt 



e'ln 



1 + Xab{0)e- 



de 



e{9) = mRe^ + $(6* - 9') In [l + Xab{0')e-'^^'^ 



dO' 



where 



Xab{e) = Ka{e - eB^)Kb{e - Ob,). 



(5.4) 



(5.5) 



It is also possible to compute the boundary entropies by using the TBA in the "L-channel" , 
as above. One finds 



In < S„|0 >< Q\Bb >= — I QabiO) In 

^TT ./_oo 



1 + e 



dO -\- est. 



(5.6) 
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where 



e(6')=mL— + / ^{6 - e')\iv I + e''^^ ^ dO', (5.7) 



and 

QabiO) = In [TZaie - eBjnb{6 - 6b,)] . (5.^ 

t CLu 



In the last expression compared to ( |4.4|) the term S{29) has disappeared because it would 
arise from LR scattering which is rapidity independent. Similarly, vanishing momentum 
corresponds now to vanishing energy and the corresponding term can be discarded. 

Despite the close similarity of the massless and massive TBA equations, what is meant 
by ultraviolet and infrared limit in the massless case is different from the massive situation. 
In fact, instead of the bulk mass parameter, we now vary the boundary energy scale 
me^^ (m just fixes the overall scale and can be considered constant) which corresponds 
physically to varying the boundary interaction. Boundary states in general will correspond 
to conformal states both in the UV and in the IR. These conformal states in general will 
have different scalar products with the conformal ground state, and the boundary entropy 
will vary along the flow. Accordingly, eq. ( f4.8|) does not hold. Physically, the boundary 
state (^]^) can now have varying scalar product with the ground state because conformal 
states can be represented as combinations of massless particle states. 

To illustrate this situation, consider the Ising model at T = Tc but with a boundary 
magnetic field. This induces a flow from the free to the fixed boundary conditions. The 
flow is described by a massless boundary matrix 

7^(^ - ^5) = -itanh ("^-^^ - (5.9) 



Using (|5.6|) and setting g =< 0\B >, one finds 

In 



In fir 



1 + e- 



dO 

— + cst, (5.10) 



3 cosh{e -In L/ Lb) 2ti 
where Lb = — e~^^. This is easily evaluated 



^ = '^V(a + l/2) U ' ^'-^'^ 

where a = L/L^. The UV limit is L — > and g — > cst-\/2, while the IR limit is L ^ cxo 
and fir — > cstl. Therefore we find g^^ /g^^ = = gf'^^^/gf'^^^'^ in agreement with 
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The above characterization of the boundary state involving massless particles is appro- 
priate for very large L. For finite L another characterization is in principle possible which 
does not involve massless particles, as we now describe. Viewing the partition function 
in the R-channel, the Hamiltonian and the Hilbert space are the same as for the confor- 
mal field theory on a cylinder. Therefore, the boundary states can be expressed in terms 
of conformal operators even though the boundary states themselves break the conformal 
invariance. Furthermore, since the conformal field theory is well defined for finite L, one 
can obtain exact partition functions for both finite L and R. For the Ising case, this was 
studied by Chatterjee 00. For perturbed minimal models, these boundary states were 
characterized in [jl^. In this section, we show how formulas of the previous sections in 
part characterize the boundary states in the simple case of the Ising model. 

Let ip{z)^il){'z) denote the left and right components of the Ising free fermion field, 
where z = t -\- ix^ ^ = t — ix. On a cylinder with x varying along the circumference of 
length L and t varying along its length R, the mode expansions read 

ip{z) = ifjr exp {—27rrz/L) 

-_ ^- _ (5-12) 

^(^) = exp {-2Tvrz/L) , 

r 

where r&'^ + l/2{r^ 2) the Neveu-Schwarz (Ramond) sector. The interpolating 
boundary condition takes the form[§]: 

^(V'-^j-iy (V' + ^)=0 (t = 0). (5.13) 
Inserting the mode expansions ( |5.12| ) into ( [5.13| ), one obtains 

where 

n - a h'^L . ^ ^. 

an = — ^ — , a = ——. (5.15) 
n + a All 

The equation (|5.14| ) implies that the boundary state must satisfy 

{^^ - a,V^_,) \B >= 0. (5.16) 



^ We obtained the results in [0] independently, in a manner presented below. 
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This fixes the boundary states up to a constants g, which can depend on L and h: 

\re'^ + l/2,r>0 J 

\B >^ = g- exp ^ a^-iA-nV^-n l^" >, 



(5.17) 



where |a > is the spin field state of dimension 1/8 in the Ramond sector. 

The factors g± depend only on L and h, and are characterized by g^ =< 0\B , 
g- =< 0\B >^. Setting the length R to zero and one of the boundary scattering matrices 
TZa to 1 in (^4.5| ), one sees that the (/-factors correspond to the boundary free energy term 



[Q term). Since the S-matrix is —1, the integral equation just implies e = mL cosh ^. Thus 
g± are just the massless limit of the expressions: 

\ogg±= hm / -z logi?(^)-7r(5(^) log(l±e-"^^^°^^^). (5.18) 

The extra minus sign in g- versus g-\- is easily obtain by including a (—1)^ in the TBA 
computation, where F is fermion number. 

One has TZ{0) = — 9), where K is given in ( |2.7| ). Defining x = mL cosh.9/2na 

and letting m — > 0, one obtains 



1 1 

I 1 , o^„.. 1 



log.,^-y^ ^log(l + e— 



[5.19) 



logg- = - / -_^log(l-e-^-«^) + -log2 

TT Jq L + X 4 



The extra term in log g- arises from the contribution 

J.™„ f'^i^O- -^("O - = i (5.20) 

(The m limit is obtained from log(l — g-'^'^coshe^ ^ log mr + log cosh ^; the log mr 
term vanishes, and the log cosh ^ term gives ( p.20| )). The above integrals are easily donell: 



2n ^a^« 21/4^2^ f«V. (5.211 



^+ r(a + l/2)Vey ' ^" " r(a + l)Ve 



The results for these integrals presented in |11] are off by a minus sign. 
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The complete boundary state is a linear combination of the NS and R boundary states: 

\Bh >= ^ {\B +sign(/i)|S >^) . (5.22) 
v2 

The partition function with magnetic fields h, h' on the boundaries is now easily computed 
as 

Zh',h=<Bh'\e-'''^\Bh>, (5.23) 
where H is the hamiltonian on the cylinder: 

i7 = ^fLo + Lo-4)- (5-24) 



L \ 12 

Using the formula ( |2.20|) , one obtains 



Zh',h = \q-'^^^9+{h)g+{h') n (l + an+i/2(/i')an+i/2(/i) 

n=0 

^ oo 

+ -s\gn{hh')g-{h)g.{h')q^'^^\{{l + a,.{h')a^{h) ^) 



(5.25) 



n=0 

where q = exp(— 47ri?/L). 

The free and fixed boundary conditions at /i = and h = ±oo respectively are con- 
formal. In these limits, the above partition function corresponds precisely to the modular 
transformation of the formula ( |2.21D . To show this, recall that in the L-channel, the 
partition function is[|5|: 

ZaB = Y.<BSixj{q}, (5.26) 
where 5" determines how the characters transform, and q = exp{—4nR/L). As h —>■ 

{h ±00), ttn = 1 (On = —1). Also, 

h^±oo: g+ = l, g. = 2+^/\ 



(5.27) 



The Virasoro characters have the infinite product expressions [|T2 



n=0 

oo 

n=0 

oo 

xi/i6(?)=?'/''n(i+^")- 



n=0 
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Let / denote free, and ± denote the fixed boundary conditions as h ±00. Then in 
these conformal limits, from ( ^.251) one recovers (|5.26D , where the only non-zero n^^ are 

1/2 1/2 1/16 1 r, . . , 

nj__|_, Uj-j^, rijjj- , Uji^, n^jj , and b is given by 

1 V2 \ 

-V2 . (5.29) 
V2 -V2 / 

(The rows and columns refer to the characters in the order 0, 1/2, 1/16.) For instance, 

Z+- = Zh'=oo,h=-oo = ^ (xo(g) + Xi/2{q} - v^Xi/iel?)) • (5.30) 
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Part II: The sine-Gordon case 



6. TBA in the L channel. 

6.1. The inhomogeneous 6-vertex model with boundary fields 

Notations closely follow those of 0. One starts from the inhomogeneous 6-vertex 
model with boundary fields h as a regularization of the boundary sine-Gordon model with 
Dirichlet boundary conditions. 

In the inhomogeneous antiferromagnetic 6-vertex model with anisotropy parameter 7, 
one gives an alternating imaginary part ±zA to the spectral parameter on alternating ver- 
tices ||T3|,p^. The scaling limit is given by taking A — > 00, TV 00, and the lattice spacing 



A ^ 0, such that R = NA remains finite. In the bulk, this provides a regularization of 
the sine-Gordon model with Lagrangian 



LgQ= / dx 



^{d(t>)^ + ^^ cos (3 SG(t) 



(6.1) 



where oc exp(— constA), = 8(7r — 7), and the field is fixed at a; = and x = R 
(Dirichlet boundary conditions) to a value that is simply related to h. 

The wave function of the inhomogeneous six-vertex model can be expressed in terms 
of a set of "roots" Oj, where j = 1 . . .n. They must be solutions of the set of equations 
1|,|16 

N [f{aj + A, 7) + f{aj - A, 7)] + 2f{a,,jH) = 

" (6.2) 

27r/j + 2^ [f {aj - am, 2'y) + f{aj + am, 2-f)], 

where Ij is an integer and all aj are positive. The function / is defined as 

/(a,6)^2ta„-'(cot5ta„ha 
and 

i7 = -/(z7,-zln(/i + cos7)). (6.3) 

7 

For h = 0, 'jH = TT — 7. By construction of the Bethe-ansatz wave function, aj > 0. Even 
though there is a solution of ( |6.2| ) with one vanishing root for any and n, we emphasize 
that aj = is not allowed because the wave function vanishes identically in this case. 
Observe that equations ( |6.2| ) are formally satisfied as well by the opposite of the roots, 
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—cxj. Often in what follows we shall consider that the roots take both signs in order to 
rewrite equations in a way which is similar to the bulk case. 



The solutions of these equations are quite intricate for arbitrary 7 [ pTTI . For simplicity, 
we restrict to the case 7 = 7 where t is an integer, and restrict to the choice e = — 1. In the 
sine-Gordon model, this falls in the repulsive regime. We make the standard assumption 
that all the solutions of interest are collections of "/c-strings" for k = 1,2. . .t — 1 and 
antistrings a ||17|| . A /c-string is a group of aj in the pattern a^^"^ — in{k — 1), a^'^^ — Z7r(/c — 
3), . . . , a*-'^^ + iTxik — 1) where a^'^-* is real. The antistring has aj = a*-"-* + ztt, where a*-"-* 
is real. 

The thermodynamic limit is obtained by sending N ^ 00. In this case, we can define 
densities of the different kinds of solutions. The number of allowed solutions of ( |6.2|) of 
type k in the interval (a, a + da) is 2N {pk{a) + p'^{a))da, where pk is the density of "filled" 
solutions (those which appear in the sum in the right-hand-side of ( |6.2| ) ) and is the 
density of "holes" (unfilled solutions). The densities pa and are defined likewise for the 
antistring. The "bare" Bethe ansatz equations follow from taking the derivative of ( |6.2|) . 
For 7 = 7r/t they can be written in the form: 

. 1 

27r(pA; + Pk) = 0^fc(") - <i>k,t-l * Pa + y^'i^kl* Pl + TTTT^fc 

1=1 (6.4) 

27r(pa + Pa) = 27r(pi_l + ptl) + ^(Wa - Ut-l) 

where * denotes convolution: 



/oo 
da'f{a - a')g{a'). 
-00 



These densities are originally defined for a > 0, but the equations allow us to define 
Pfc(— a) = Pk{c() in order to rewrite the integrals to go from —00 to 00. The kernels in 
these equations are defined most easily in terms of their Fourier transforms 

/oo 7 I 1*00 

— e^"*"/"/(a), /(a) = - / e-^"*"/"/(x)dx. (6.5) 
27r TT ./_^ 



One has 



7 , X e ^coshxsinh(t — /c)xsinh/x ,^ ^, 

^kAx)=dah-2 T-r r-r- , (6.6) 

smh X smh tx 
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for k > I with (pik = (pkh and 



^ smh(t — k)x 

ttk = cos Atx/TT 

sinn tx 

^ ^smh{t — H)x smhkx sinh(t — 2/c)x/2 ^ 

sinhxsinhtx sinhte/2 (6-''') 

^ ^sinhi/x sinh(t — 2)x/2 ^ 
" sinhtx sinhtx/2 ' 

with in particular ai{a) = \ /(a + A, 7) + /(a — A, 7) ,(j)ii{a) = —/(a, 27). The 
boundary manifests itself in the first term in Uk] notice that even for /i = 0, it still 
modifies the equations. A few technicalities account for the other terms (these are relevant 
here because we are interested in subleading boundary effects). The second term in Uk 
arises from the fact that the sum in ( |6.2| ) does not include the term m = j; the integration 
over densities includes such a contribution and so it must be subtracted off by hand. The 
third term in Uk arises because p and p'^ are defined for allowed solutions, while as already 
explained, a = is not allowed because it does not give a valid wavefunction. Since it is 
a valid solution of (|6.2| ) but is not included in the densities, we must subtract an explicit 



^d{a) (corresponding to 1/2A^ in Fourier space). Explicitly, one has 



ui = 2/(ct, 7if) + 2/(2a, 27) - 2Tvd{a). (6.8) 



For compactness we rewrite (|6.4| ) as 



2TTe^''\pk+pt)=ak{a) + J2<Pki*Pi + ^Uk, (6.9) 



2N 
I 



where e*-'^-' = — 1 for the antistring. 

The energy reads, with proper hamiltonian normalization. 



E 



latt 



2N 



1 f°° 

/ ak{a)pk{a)da. (6.10) 



It is easy to write the thermodynamic Bethe ansatz for this model. One finds that 
the TBA equations, since they are obtained by a variational method, do not depend on 
boundary terms, and read as usual 

A, 

I 
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-afc(a) =Tln(l + e^^)-TVe(')-^*ln(l + e-^0: (6-11) 



where 

Aki{a) = 2'Kt^''Hu5{a)-^ki. (6.12) 
The free energy does depend on the boundary term and reads 

k 

rji /»00 

^ J-oo 



-ilatt 



da 
2^' 



(6.13) 



In the above formulas the temperature T corresponds in the two-dimensional point of view 
to having a cylinder of radius L = 1/T. We can deduce from this result the ground state 
energy. Indeed recall that the ground state is obtained by pfc = 0, /c 7^ 1 and = so 

j^iatt ^ _^ r ^^(^)|^-|^ - ^ n,\e^\P, (6.14) 



00 27r 2 2tx 



where from ( |6.11|) we have 



er = --'"'^^f/"^ (6.15) 
t cosh X 



Replacing and using (|6.7| ) we find 



Eiti = - - / COS 



^ 2 ( sinh(t-l)x 



^"''^ TX j _ao V TT / coshxsinhtx 



j^iatt r cos(Ate/7r) /^ sinh(t-g)a: ^ sinh(t - 2)a;/2 

^'^^y 2n J-oo cosh a; \ sinhtx sinhtx/2 

where we used the formula 



(6.16) 



/oo noo 
a{a)b{a)da = 2t a{x)b{—x)dx . 
-00 J —00 



6.2. Boundary energy and entropy of the SG model with Dirichlet boundary conditions 

In the continuum limit A — 00 the energy contains various terms. We keep only the 
finite part which is obtained by closing the above integrals in the upper half plane and 
selecting the pole aX x = i^, leading to 

Ebuik 'Y 

m ( sin(t-if)7r/2 tyr ^ ^^'^^^ 



Ebdry = - TT 2 7^ A 

2 \ smtn/2 4 
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where m is the sohton mass, 

m = 2e-*^/^ (6.18) 

All these results trivially generalize to the case of two different boundary fields by splitting 
the H dependent terms into the sum of an H and an H' term. The bulk result agrees with 
what is obtained by other methods For completeness it might be useful to write 



the boundary energy in terms of the parameter ^ of 

^ m /cosft — l)f cosft — l)f ' .^ 

Ebdry = -IT .\ ,^ + , - cott7r/4 - 1 . 6.19 

2 V smt7r/2 smt7r/2 J 

As in the bulk case, when t is even, there are additional logarithmic terms. 

We can as well take the continuum limit of the lattice Bethe equations to get the 



TBA equations for the sine-Gordon model. This is explained in details in and ||2C1|| . 
We simply notice here that in the IR limit, all e's go to infinity so one finds a vanishing 
boundary entropy in ( |6.13| ). In the UV limit, it is easy to solve the system (|6.11| ). One 
finds 

' Hk + 2) (6.20) 

1 H~ Xa 

where Xk = e~^^. Now we have the two identities just proved by inspection 



k=i 

t-2 



k{k + 2) t 
(A; + 1^2 

K ill 

k = l 

Using the fact that 



(6.21) 



u,{0) = k[2---—j 



one finds that 



^I,_i(0) = (t-1)(2-^-^) , (6.22) 
^ . X 2 2H 

Ua{0) = -2+- + —, 

/oo 
Uk{a) = 0. (6.23) 

Hence the boundary entropy is actually the same in the UV and IR limits, as was argued 
for a scalar theory ( [1.1| ). 
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7. The Destri De Vega equations for the boundary sine-Gordon model 
7.1. The DDV equations with boundary conditions 

We now would like to compute the complete Casimir effect in a theory with boundary. 
TBA in the R channel is pretty intricate in the non diagonal case. We use an alternative 
method elaborated by Destri and De Vega (DDV) in the periodic case. The lattice regu- 
larization plays a crucial role here - it is actually not known yet how the DDV method is 
related to the standard description of the theory using excitations and S matrices. Inter- 
estingly, the result are very close to those of the TBA of scalar theories in the R channel, 
and with some more work could probably be considered as a more rigorous proof of the 
results presented in the first part. 

For that consider eq. (^]^) which we rewrite as 

2Np{aj) +Pbdryiaj) + ~" ""^^ ^ ^^^J' ^'^'^^ 

am>0 

where the sum runs over all roots (including m = j) and we introduced the notations 

= ^ [/(« + A,7) + /(«- A,7)], PH{a) = f{a,-fH), (f){a) = Ma) , (7.2) 

Pbdryia) = 2pHia) - (j){2a). 

The ground state is obtained by filling the real positive solutions, aj = excepted. 
This corresponds to the choice Uj = 1,2, . . .. Recall that if aj is solution of ( [7.1| ) with 
some Uj, so is formally — (with —Uj). Given the set of roots {aj > 0} representing the 
ground state, one can construct the counting function as follows: 

f{a) = 2iNp{a) + ipbdryioi) + i ^ (/>(«- a^) +(/>(« + a^), (7.3) 

a77i>0 

Define then 

y(a)=e-^("\ (7.4) 

We have Y{aj) = Y{—aj) = 1 for every root aj from the ground state, as well as y(0) = 1. 
Therefore we can rewrite ( [7.3|) as 

f , Yia') da' 

f{a) = 2iNp{a) +ipbdryia) -i(p{a) - (p{a-a')- tft^t;-, (7.5) 

In 1 - Y (a'] Zn 
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Figure 3. Contours of integration in the DDV method. 



where at the right-hand side takes care of the unwanted contribution of the pole 

a' = and the contour C consists of two parts as shown in figure 3, Ci above and C2 
below the real axis. 

Like in the bulk case 0, simple manipulations allow us to rewrite this non-linear 
integral equation as 



f{a) = 2iNP{a) + iPbdry{a) + [ $(« - a') In - e^^"')) da' 
+ f ^{a - a') In (l-e-f^'"' A da'. 

J Co 



'C2 

In equation ( |7.6D one has 



2tx^ J _^ 2 sinh(t — l)x CO 

together with 



cosh X ' 



(7.6) 



(7.7) 



^, , ;+°° , e-itea/TT _ I cos(Ate/7r) , 

Pia) = / dx ' 7-8 

3 —IX 2 cosh X 
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and 

aa; ■ 



bdry 



-IX 



smh(t — l)x cosh a; 

(7.9) 

sinh(t — 2)x/2 coshtx/2 smh(t — 2)x 

sinh(t — l)x cosh x 2 sinh(t — l)x cosh x ' 

Obtaining ( [7 .61) requires some care with the definition of logarithms. One proceeds as 
follows. Before integrating by parts in the integral over C2 one factors out ^ \n.Y{x): 

^^""^ in[i-y-i] + -^iny(x). 



1 — Y{x) dx dx 
Then both integrals over Ci and C2 can be taken by parts, resulting in 

/■+°° • da' 
f{a) - I (j){a - a')f{a')— = 2iNp{a) + ipbdry{a) - i<f>{a) 

J-00 

/'+°° • da' 

- / (j){a-a' -iO)ln[l-Y{a' + iO)] 

J-00 27r 

/+00 rlru' 
^{a -a' + iO) ln[l - Y'^a' - iO)] , 
-00 27r 

(surface terms from Ci and C2 cancel against each other provided A^, A are finite) . To 
make the source term vanish at infinity we take the derivative of both sides in the latter 
equation, after which it can be Fourier-transformed and "dressed" by the factor (1 — ^)~^. 
Finally, one goes back in the rapidity space and integrates the equation, using /(O) = 0, 
to obtain ( [7.6|) . 

The energy of the ground state configuration can be expressed as 



2 

j^latt 

1 • , . , 1 /■ A . Y{a) da 



+ l [ Ma- A)- ... , 

t t Jc 1 ^('^) 2«7r 

where /(a, 7) = f-yia). One finds after exactly the same manipulations as above, 

1 • 1 /■+°° da 

E = -MA) - - / /;'(a - A + zO) ln[l - y(a + zO)] — 

»+oo 1 r+00 



1 r^"^ . da 1 r^°^ 

+ - / f'^{a-A- zO) ln[l - Y-\a - .0)]— + - / f'^{a - A)f{a) 



(7.11) 



2m' 
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Substituting (|7.6|) instead of f{a) in the last term of the latter, we obtain 

j^iatt ^j^iau ^ j^iatt _ \ f^^ siy - A + ^0) ln[l - Y{y + tO)] ^ 



(7.12) 



+ -J_ s{y-A- zO) ln[l - Y-\y - ^0)]^, 
where we defined s{y) by 

2cosh7A; " 4 cosh(ty/2) " ^'^ ^^ 

The last two terms in (|7.12|) represent finite-size corrections to the ground state energy, 
while 



'bulk -^bdry —jJ-yV^) ~ J / ^6 J-y 



- ^" 27r-(/>(A;) (7.14) 

1 ■ 2A^ /'+°° ■ 

= -7f'yi^) 1- f^{a- K)p{a)da, 

^ ^ J —00 



where the function p(a) defined so satisfies the following equation: 

27rp(a) = p{a) + </> * p{a) + ^ (]36dry(a) - </>(«)) , (7.15) 

which can be checked by solving this linear equation in Fourier space. Introduce pi(a) = 
p{a) - 5{a)/2N. Then 



2N 
1 



+00 



Eluik + Eity = / /7(« - ^)pi{cx)da (7.16) 



00 



and, by virtue of (|7.15D, pi satisfies the equation 



27rpi(a) = p{a) + ^ * Pi(a) + ^ {Pbdry{oi) - 2n6{a)) . (7.17) 



'^latt 



Hence pi is the density of the ground state configuration ( p^ ) (see also (|0|)) and -E^^/^ 
and -E^^^y coincide with the quantities computed in the previous section. 
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1.2. The continuum DDV equations 

Having checked the correct values of bulk and boundary energies, we now let the 
cutoff A ^ oo according to A = | log(2/mA) with the size of the system R = NA and the 
physical mass m = -^e"*^/^ fixed, and work only with the renormalized theory. In that 
limit one has, evaluating all integrals in Fourier transforms and keeping the leading terms, 

t^ 

s(a + A) + A) = -— msinh^, iVP(^) = mi^sinh^, (7.18) 

where we set 

Recall indeed that when one studies the excitations of the model, a relativistic dispersion 
relation is obtained provided the rapidity in the relativistic theory and the "bare rapidity" 
of the Bethe excitations are related by ( [7.19|) . We redefine implicitely all functions to 
depend on 9 from now on. The energy therefore reads now. 



E =Ebuik + Ebdry + ^ f m sinh e\n(l - e^^^^^ ^ 
+ rnsinh6'ln (^1 - e"^*^^) 



where / is solution of the integral equation 



(7.20) 



f{e) =2imRsmhe + iPbdry{0) + J $(6* - 6') In (^1 - e^^^'^) dO' 



(7.21) 



and 



27r2cosha;sinh(t-l)x ' ^ ' 



which can be identified with 



$(0) = --^-|ln^++(^), (7.23) 
2z7r du 



where S-\-\- is the soliton-soliton S matrix element. We have also from (|7.9|) and the identity 
(18) in [0 (see also J2T[) 

PbdryiO) = -- lnS++{2e) + - ln7^++(^) + n, (7.24) 

(so -jgPbdry = + 27r5(6')). The additional tt in the above equation arises because the 
has an overall minus sign [|22| (eg 5" = — 1 in the free fermion case t = 2). 
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1.3. The Casimir effect 

We define the effective central charge by the formula 

^Ceff 



E = Ebulk + Ebdry 



24R 



(7.25) 



Observe that in the approach of the first part, we dealt with a renormalized theory with 



the vacuum energy subtracted. The last term in (|7.25|) is a renormalized energy, and it is 
what we will compare later to the results of the first part. From ( [7.20| ) we find 

6R 



Ceff{mR) 



m sinh ^ In ( 1 — e 



Ci 



d9 



+ J msinh^ln (^1 -e--''(^))| 



d9. 
(7.26) 



To study the ultraviolet behaviour of the above expression, let us use f{6) = f{6) and 
rewrite (|7.21| ) as 



f{0) = 2imRsmhe + iPbdry{0) - 2i 
and ( |7.26| ) as 



^{e~e') Imln 



CeffimR) 



12mR 



d9 sinh 9 Imln 



1 



1 _ efis'+io) 



J(9+iO) 



(7.27) 



(7.28) 



It might be useful to remind the reader of the form of the corresponding bulk equations 

"l + e/(e+.o)- 
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[mR) 



6mR 



with / satisfying 

f{9) = imRsmhO + iuj — 2i 



sinh^ Imln 



^{6-6') Imln 



1 + e 



d9, 



f{e'+io) 



de' 



(7.29) 



(7.30) 



and uj is the twist of the 6-vertex model. Note the deep similarity between these two 
systems; the factors of 2 can obviously be absorbed in a redefinition of mR and the minus 
sign in the arguments of logarithms in a redefinition of /, so the only essential difference 
is that the twist angle u is replaced by Pbdry It is well known that the twist corresponds 
to a soliton fugacity [0 so we recover the result of the first part that the boundary acts by 
some effective, rapidity dependent fugacity. 

In the limit when i? ^ 0, only the region \6\ large contributes to Ce//. Let us focus 
on the limit ^ >> 1, the results for negative 6 following by symmetry. Then one finds 
f ^ Ik, where 



fK{0) = imRe^ + iPbdry{oo) 



2i 



^{9-9') Imln 



1 _ e/i^(e'+io) 



d9' 



(7.31) 
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together with 



6 

CKimR) = I mRe^ Imln 



1 _ g/x(e+iO) 



de. 



(7.32) 



It is now useful to recaU some weU known results about dilogarithms |2J]. Define 



Lix) 



du 



ln(l + w) In w 



u 



(7.33) 



Assume 



/oo 
dy G{x- y) Im In [1 + ^(y + iO)] , (7.34) 
-oo 



with G an even function. Then one has 



Im / dx(t)'{x) In [1 + F{x + ze)] = -Re {L[F(-oo)] - L[F(oo)]} 

J— oo 2 

+ ^Im {(/)(oo) In [1 + F{oo)] - (/)(-oo) In [1 + F(-oo)]} , 



(7.35) 



(where we did not write the iO part of some arguments for simplicity). Set F = e-^^~^'^ 
and denote Pbdryi^) = cr. Then, according to ( |7.31J ): (p = mRe^ + a — tt, G = — We 
have (j){—oo) = a — tv, 4>{oc) = oo. One has also F(oo + iO) = 0, and from this and ( |7.31j ) 
one can get the value of F at — oo: 



/k(-oo) 



=0" — 2 Imln 

t - 2 



I — g/x(-cx)+jO) 



^{e)de 



--a + 



t - 1 



Im In 



I — QfK{-oo+iO) 



(7.36) 



So, if e^^(-°°+^o) = e^'^ one may use Imln(l ± e'^) = j: ln(±e^^) to find 



t — 1 TV 

exp <^ 2i — - — a + 2zy ^ . 



(7.37) 



From ( |7.35| ) it follows that 



Im / mRe^ ln(l - e^^^^+^o))^^, _ 1 ReL(-e^"') - - (a - tt) ("^ — -a - tt + - 



In the region ^ << 1 we have /i ~ /a? and similar calculations yield 



' 1 1 

Im / mi?e-^ ln(l - e^^(^+*°))d^ = — ReL{-e-'^) + - (tt - a) 
-oo 2 2 



t-1 n 

a + 71 

t t 
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Collecting both ^ >> 1 and ^ << 1 contributions we obtain 

"7r2 t-1 



6 r 1 



L{-e''^)) + L{-e-'''^) +- — - {a - nf 

t 



6 



6 t 



{a - ttY 



From ( |7.9| ) we get 



Finally, from this we find 



Pbdry (OO) = a = 271 



nH + H' 
2 t-1 ' 



t-1 f H + H' 
cuv = 1 - 6 — ■ — I 1 



t V 2(t-l) 

In the case with no boundary field, H = H' = t — 1 so cuv = 1 as expected. 



(7.38) 



(7.39) 



(7.40) 



Remark. So far we tacitly assumed that < H < t — 1. In general, from relation ( |6.3| ) 
follows that when /i > 0, H varies between —t — 1 and —1, while when /i < 0, — 1 < 
H < t — 1. To generalize our results, we should use the most general form of pn'- 

PH{k)= r pHia)e'''"da = 2n sign{H f''^^^!'~ "^^"^^ , -nK'jHKn, (7.41) 
J-oo smhnk 

where we defined ujh = It-^I- For —2tv < •yH < —ir set uh = 2% + •yH and sign(ii/') = 1 
in ( [TilD . Then (|7]3|) generalizes to 

2 TT — 7 

if H and if' are both positive or —2t < {H, H') < — t, and 

A-K - Uh - i^'n 

(J = ZTT 

2 TT — 7 

if they are both negative, but greater than —t. In the case when ^ < H < t — 1 and 
-t < H' < -1 (that is, h<0,h' >0) we get: 

TT - - 27 

a = TV -\ . 

2 7r — 7 

So, when — = 27 we have a = tt and cuv = 1, as in the free case. The condition 
iv'jj — Uh = 27 is equivalent to h = —h', as could be seen from ( |6.3| ). That c = 1 when the 



two surface field are real and opposite is well known from lattice studies pH . 
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8. Remarks on the DDV equations for minimal models and excited states 



A particularly interesting case is when the XXZ chain or the inhomogeneous 6-vertex 
model commutes with the quantum group Uqsl{2). In that case h = —h' = 2zsin7 and 
the net result is that all H dependent terms simply disappear from the equations so in 
particular 

m, / tir \ , ^ 

8.1 



bdry 



m 

~2 



tn ^ 
cot — + 1 
4 



At the N = 2 supersymmetric point, t = 3 so the boundary energy vanishes, a result well 
expected from supersymmetric considerations. More generally, it vanishes if t = 4n + 3, n 
an integer. Notice that the bulk energy vanishes for t an odd number (as a consequence 
of the generalized fractional N = 2 supersymmetries studied in [|26|). 

In the quantum group symmetric case one has H + H' = 2t so from ( [7.40|) 

6 



c = 1 



t{t-l] 



(8.2) 



the expected result for the restricted sine-Gordon model P] [p8| [p9|. 

In the quantum group case, all -ff-dependent terms simply disappear from the equa- 
tions (in notations of this corresponds to ^ — > ±ioo). We find then the value of the 
boundary matrix element 



ln7^ 



dx sinh(3x/2)sinh(t-2)x/2 . 2 ^ 

n ^ — ; n ^ sm —xO . 

X sinh(t — l)x/2 sinh(t — l)x tt 



(8.3) 



In the last equation we suppressed the label ++ because the scattering of solitons and 
antisolitons can be made identical by the same change of gauge which ensures commutation 
of the bulk S-matrix with Uqsl{2). To understand ( |8.3D , we can use the reformulation of 
the quantum group symmetric bulk-scattering in terms of the 0{n) model. In that case, 
the bulk S-matrix can be rewritten ||3^ 



S. 



p{-ie) 



ab 



[sinh A(z7r - 9) S^J^ + sinh X9 daj^'^] 



(8.4) 



where 



p(—i9) = —i sinh ^ Xiin — 9) exp 



dx smh ((y - l) x) 2 ^ 

r^T -TT sm -x9 

X sinh (I j cosh (x) tt 



(8.5) 



The labels correspond to colors in an 0{n) symmetric model, running formally a = 1, . . . , n. 
One has n = 2coshz7rA so the sine-Gordon model corresponds to n < 2 and one must 
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continue things analytically in n. In our previous notations, A = In this approach, 
the particles are self conjugate so the unitarity crossing relation reads 

K - o) = s'J:,{2e)nl {i^ + e) . (8.6) 



We deal with the 0{n) model with free boundary conditions so by 0{n) symmetry, TZ 
actually does not depend on a let us just call it TZ. Then ( |8.6| ) reads 



(8.7) 



where we used n = 2 coshzTrA. The latter equation appear in [^, and its minimal solution is 
precisely ( |8.3| ). More results about the boundary 0{n) model will be presented elsewhere, 
together with a kink interpretation for minimal models. 

Now in [^, it was shown how one can obtain the various conformal weights hrs 



n(^i^-9^ = pLJ^ [n smh2Xe + sinh X{in- 29)] n[i^+ 9 
= sinA(7r - 2i9)p{-2i9)n (^i^ + 9 



by changing the boundary conditions and the total value of the spin in the XXZ chain. 
The change of boundary conditions proposed in was to project the first r spins near 
a boundary onto the spin r representation of Uqsl{2), which leads to characters hr,i. 
Now let us recall (see for discussion, see also ||3^) that if 7^_|__|_, 7^__ are boundary 
Refiection-matrices satisfying unitarity and crossing, one can obtain new solutions of the 
same equations by defining 

nf+{9) = a{9 - 9o)a{9 + ^o)7^++(^) 

7^L^1 {9) = b{9 - 9o)bi9 + 9o)n--{9) + c{9 - 9q)c{9 + 9o)n++i9), 

where a, 6, c denote the standard elements of the bulk sine-Gordon S-matrix, and 9o is 
a free parameter. The process can be iterated again to give R-matrices TZ^'^\ Now one 
can again let ^ ±zoo to obtain what we conjecture to be the continuum limit of the 
projection discussed in [p5|. The DDV equations read as above, but with TZ^'^^ replacing 



TZ. The effective central charge follows from P^dry {oc). Using 

lnS'_|__|_(cx3) = —n- 



i ' t-1 

one finds 

Pu^^.Jno) = TV ( 

t-1 ' 



Pbdryioo) = TV ( ^-^ - (r - 1)^-^ ) . (8.9) 
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Using the equations (|7.36|) to ( p.2| ) this leads to 



t - 1 



t 



) 



2 




1-6 



t 



1 



t - 1 



r 



c - 2Ah. 



rl- 



(8.10) 



Hence we believe we can observe the conformal weights hri using these "boundary fused" 
R-matrices (what value to give to 6q requires more discussion). 

Changing the value of s is more difficult from the DDV point of view. On the lattice, all 
one has to do is to look at non vanishing spin m. Then one observes /ir,i+2m5 corresponding 
to 



Unfortunately all our equations are defined modulo 27r. So it is possible to observe /ir,2 
only. For that purpose simply subtract tt to (|7.24| )(this corresponds to having a lattice 
with an odd number of sites). 

9. The massless DDV equations 

The structure of the DDV equations is pretty transparent in continuum terms: the 
function $ takes care of the interaction between solitons and antisolitons, and the term 
Pbdry takes care of all the boundary effects. It is then natural to guess what the DDV 
equations would look like in other cases. A case of special interest is the bulk massless 
case. If the boundary is also conformal invariant, one gets the set of equations ( [7.31| ) and 
( [7.32| ), and we see that the particular conformal boundary condition is fully characterized 
by a rapidity independent phase. If the boundary is perturbed and the perturbation is 
integrable, its action can be described by a boundary scattering. If moreover this scattering 
is diagonal and identical for solitons and antisolitons, it is easy to write the equations as 



Here 9b characterizes the strength of the boundary perturbation, and me is the usual 
"Kondo" or "boundary" temperature. The only 6 dependent term in Pbdry is the one 
involving the boundary S-matrix. The simplest case corresponds to the following boundary 
S-matrix 




(8.11) 




(9.1) 




(9.2) 
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corresponding to an interaction of anisotropic Kondo type on one side of the strip, and we 
chose the usual Dirichlet boundary conditions on the other side. Then one has 

PbdryiO -dB) = ] Intanh - . (9.3) 

As a result, the energy scales as the ground state of the free boson theory with c = 1 in the 
UV, but according to (|7.40|) in the IR one observes an excited state of conformal weight 



t-1 f Ptdryioo)\'^ t-1 



At \ n J At ^ ^ 

Microscopically, the IR limit corresponds to having an infinite boundary field in the XXZ 
chain and one easily checks directly the validity of ( |9.4|) . 



10. Conclusions 



To conclude, consider the energy of the renormalized theory, E 



24i? 



rf. So far, the 



relation between the DDV equations and the underlying scattering structure is unknown. 
Nevertheless, it is likely that the two terms (Ci and C2) in the above equations correspond 
to the existence of "soliton" and "antisoliton" . By naive extension, the result for a scalar 
theory should have the same form but with only one of these terms. After shifting the 
countour Ci to Im(^) = —in/2 and changing a few notations, one is thus led to the 
following conjecture for the ground state energy of a scalar theory with bulk scattering 
given by $: 



E 



— — / cosh 9 In 
Att 



1 + e" 



d9. 



(10.1) 



where now 



e{e) = 2mRcosh9-iPbdry{0 + iT^/'^) + iT^+ / ^{9 ~ 9')Im\n 
and 



l + e-^(^') 



m 



1 d 

2hTd9 



lnS{9). 



d9', (10.2) 



(10.3) 



Observe that, introducing the notation K{9) = TZ — 9) we can rewrite 



iPbdry{9 + in/2) = -lnK{-9)K'{-9) -liiS{29), 



;io.4) 
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where we used the relations S{9)S{—9) = 1 and S{iTT — 9) = S{9). We can also use the 
boundary crossing relation 0] 



K{9) = S{29)K{-9), 



(10.5) 



together with K{9) = K{—9) to rewrite 



-iPbdry{9 + iTv/2) + in 



\nK{9) -\nK'{9). 



(10.6) 



This is exactly the equation we deduced in section 3. 
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